Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



Absorption and scattering by interstellar dust: 
XMM- Newton observation of Cyg X-2 

E. Costantini 1,2 , M. J. Freyberg 3 & P. Predehl 3 

1 SRON National Institute for Space Research, Sorbonnelaan, 2, 3584CA, Utrecht, The Netherlands 

2 Astronomical Institute, Utrecht University, P.O. Box 80000, 3508TA Utrecht, The Netherlands 

3 Max-Planck-Institut fur extraterrestrische Physik, Giessenbachstr. 1, D-85748 Garching bei Miinchen, Germany 

Received / Accepted 

Abstract. We present results of the XMM-A r euiton observation on the bright X-ray binary Cyg X-2. In our analysis 
we focus upon the absorption and scattering of the X-ray emission by interstellar dust distributed along the line 
of sight. The scattering halo around Cyg X-2, observed with the CCD detector EPIC-pn, is well detected up to 
~7 arcmin and contributes ~5-7% to the total source emission at 1 keV, depending on the dust size distribution 
model considered. For the first time spatially resolved spectroscopy of a scattering halo is performed. In the halo 
spectrum we clearly detect the signature of the interstellar dust elements: O, Mg, and Si. In the 0.4 — 2keV 
band, the spectral modeling of the halo shows a major contribution of silicates (olivine and pyroxene) . The spatial 
analysis of the halo surface brightness profile shows that the dust is smoothly distributed toward Cyg X-2 at least 
for ~60% of the path to the source. However, given the substantial pile-up, we could not investigate fainter or 
narrower components of the halo. Within this observation limits, the data do not show preference for a specific dust 
size distribution; namely the Mathis, Rumpl & Nordsieck (1977) or the Weingartner & Draine (2001) model. In 
this analysis we used the Mie theory to compute the differential scattering cross section. The RGS data were used 
to investigate the ISM absorption. In particular the absorption spectrum shows complexity around the oxygen 
edge at ~0.54keV, which cannot be explained in a unique way: absorption by molecular oxygen or ionized atomic 
oxygen, as proposed in other studies of Cyg X-2. Combining the RGS results with the additional information on 
dust grains provided by the EPIC-pn spectrum of the scattered radiation we estimate a column density for dust 
absorption by oxygen, provided that it is locked in silicate grains. 
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1. Introduction 

The observed light from a source is obscured by the in- 
terstellar matter (ISM) through the combination of two 
processes: absorption and scattering. Absorption is due to 
both gas and dust, whereas scattering is attributed to dust 
alone. Differently from the IR to UV wavelengths range, in 
the X-ray regime the observation of absorption and scat- 
tering by interstellar dust (ID) are strongly coupled. Thus, 
in the X-ray regime, the two extinction mechanisms can 
be simultaneously observed and studied. If an X-ray emit- 
ter is located behind a layer of dust, its radiation will be 
absorbed and at the same time scattered into the direction 
of the observer. In the X-rays the scattering mechanism 
is no longer explained by the simple Rayleigh formula. In 
particular the scattering angle is in this case very small 
(#scatt fx (V a ) ^ 1°)' forward directed, dependent on the 
wavelength of the incident photon (A) and the size (a) of 
the grain. The small scattering ang le results in a h alo of 
diffuse emission around the source <|Overbecdll965|h The 



energy range in which absorption and scattering can be 
studied, is a strong function of the equivalent hydrogen 
column density of the medium (Nn). Indeed the X-ray ra- 
diation is obscured by absorption depending on the value 
of JVh: I = lQe~ NrL<T , where Iq is the source radiation and a 
is the absorption cross section. Through the analysis of the 
absorbed spectrum, information on the chemistry, column 
density and abundances of the ID grains can be inferred. 
Simultaneously, the spectral and spatial properties of the 
X-ray halos can be analyzed. The halo intensity, angular 
extension, and spectral distribution are a function of the 
size distribution and composition of the scatterers (the 
dust grains), their distribution along the line of sight, and 
the spectral properties of the source illuminating them. 
Sources with faint halos have a hydrogen column density 
which is low enough to not completely absorb the soft X 
rays. High sensitivity instruments are needed to study the 
emission of the scattering halo, which is very weak com- 
pared to the brightne ss of a background sourc e (up to 20% 
of the soft emission, IPredehl fc Klose1ll996j) . Faint halos 
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are important to study ID chemistry. Indeed the energy 
range where scattering occurs (approximately 0.3 — 2keV), 
includes possible features of ID components, primarily 
oxyge n (0.54 keV), magnesium (1.3 keVV, and silicon (1.84 
keV) l|Predehl fc Klose Ill99rl iDraine 1120031 hereinafter 
D03). An energy resolution of ~ 80 — 150 eV over the en- 
ergy band of interest for the scattering process, makes it 
possible to spectroscopically investigate the features of a 
faint halo. Previous studies of scat tering halo profiles were 
carried out with Einstein (e.g. iMauche fc Gorenstein I 
Ij983 lOallagher Gash fc Creenl ll 995h and ROSAT 
llSmith fc Dweklll998t Ipredehl fc Schmitt I Il995i here- 
inafter PS95). According to those findings, the ID 
size distribution appeared t o be consistent with the 
iMathis. Rumpl. fc Nordsieck I l|l977l) (hereinafter MRN) 
model. The MRN model includes a mixture of car- 
bonaceous and silicate materials, with size distribu- 
tion a~ 3 - 5 , for 0.001 < a < 0.25/zm. Alternative 
grain size distribution models have been proposed, dif- 
fering mostly in the chosen boundaries of the grain 
size range, the slope of the distri bution and the in- 



ner structup3^of_&egraiii itself (Mathis & Whiffcn 
Il989t lLandgraf et alJl2000t IWeingartner fc Draine Il200ll 
hereinaft er WD01). From R OSAT data (Nova V1974 
Cygni, IDraine fc Tan I 2003) and the Chandra obser- 
vation of GX 13+1 l|Smith. Edgar, fc ShafeTI 120021) the 
role of grains of size much larger than 0.4/Ltm has been 
stated not to play a major role, at least in the diffuse 
ISM, contrary to the Solar System environment, where 
grain s up to 1 um should have a significant contribu- 
tion <|Witt. Smith, fc Dwek ll200lL On the other side of 
the range of the grain size, a better understanding of 
the infra-red spectrum of ID allowed to state the im- 
portance of ultra-small particles, the Aromatic Polycyclic 
Hydrocarbons (PAH), with radius a < 30 A. Such a con- 
tribution was included in the WD01 model. 

Finally, it has been recognized that a simple analytical 
computation of the differential scattering cross section, the 
so called Rayleigh-Gans (hereinafter RG) approximation 
, could be misleading if applied to halo energies < 1 keV 
and / or large grains a > 0.25/im. The full Mie theory IjMie I 
Il908h . from which the RG approx i mation is derived, ha d 
to be used l|Smith fc Dwek 111 1 9981 IDraine fc Ta"n~ll2003h . 
The models applied to X-ray scattering halos are the re- 
sult of a deeper knowledge of ID properties gathered at 
longer wavelengths. Due to the low resolution of early X- 
ray instruments, only integral properties of the dust could 
be studied, adding relatively little information on the na- 
ture of ID. On the contrary, with the X-ray observatories 
now flying, we can address other issues like: (i) the chem- 
ical properties of dust particles that scatter X-rays, (ii) 
abundances and depletion in the ISM, and (Hi) the actual 
distribution of dust along the line of sight. 
In this paper we present the RGS and EPIC-pn analysis of 
Cyg X-2, located at Galactic coordinates I — 87.33°, b = 
— 11.32°, behind a dust layer with equivalent TVh column 
density of the order of ~ 2 x 10 21 cm~ 2 , which produces 
a relatively weak scattering halo. This makes Cyg X-2 an 



ideal candidate to study both the spatial and spectral dis- 
tribution of the halo at energies softer than 2 keV. The 
fractional halo intensity of Cyg X-2, defined as the inten- 
sity of the halo extended emission over the total observed 
emission, was estimated from ROSAT-PSPC to be 3.9% at 
1.06 keV (PS95). Now, with the high sensitivity of XMM- 
Newton, the halo can be resolved and analyzed down to 
0.4 keV with the EPIC-pn. Absorption by the ISM toward 
the line of sight of Cyg X-2 was studied with the R GS. 
Recently Cyg X-2 was stu died bv iTakei et alJ (2003), us- 
ing Ch andra-LiETG and bv l.Iuett. Schulz. fc C hakrab artvl 
(2004) using Chandra-WETG. In each analysis, the absorp- 
tion features in the spectral re gion of the oxygen e dge were 
interpreted in different wavs. Tlakei et all l)2003l) claimed 
to have detected absorption by mole cular oxygen, while 
l.luett. Schulz. fc Chakrabartvl l)2004l) interpret those fea- 
tures in terms of mildly ionized oxygen in the ISM. 

The paper is organized as follows: In Sect. 2 the princi- 
ples of the scattering halo theory are presented. In Sect. 3 
the analysis of RGS and EPIC-pn data of Cyg X-2 is 
shown. Sect. 4 describes the careful extraction of the infor- 
mation on scattered radiation. Sect. 5 is then devoted to 
the spatial and the spectral modeling of the scattered halo. 
Finally, in Sect. 6 we discuss our results, and in Sect. 7 
the conclusions of this work are shown. 



2. The Halo Theory 

The intensity of the light scattered by dust at a scattering 
angle # sca , assuming spherical grains and single scatterin g 
has the following general form (e.g., IMathis fc Lee Il99llh 



1(0; 



„, a , = / F(E)dE / n(a)da x 



r 1 f(x) da 

, (i~x) 2 dn 



dx. 



(1) 



F(E) is the spectral energy distribution of the 
source; a is the grain radius with number density n(a); 
x is the fractional distance of the total path, from the 
source (x = 1) to the observer (x = 0), at which the 
scattering occurs; and f(x) is the normalized spatial 
distribution of the scattering sites. Since we consider dust 
that is evenly distributed, f(x) = 1. For the scattering 
angle 6* sca , it holds that # sca = o bs/(l — x )- The term 
da/dQ. is the differential scattering cross section, a 
function of 8, a and E. There are two critical terms in 
Eq ^ One is n(a), which depends on the physical and 
chemical state of the dust grains. The other crucial term 
in evaluating the scattered emission is the differential 
scattering cross section for which the exact solution 
is giv en by the Mie theory (|Mie I Il908t Ivan de Hulst I 
119571) . This describes the scattering and absorption of an 
electromagnetic wave by spherical solid particles. 
Quantitatively, the refraction index m of a given material 
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can be written as (e.g., iHenke. Gullikson. fe Davis 1993): 



ottering Cross Section 



r A 

m=l- -^-^2n q f q (0) 



(2) 



where r c is the classical electron radius, n q the num- 
ber density of the atom q. The quantity f q {0) is the 
complex forward atomic scattering factor for the ele- 
ment q. This definition only holds for small scattering 
angles, as is the case in the X-ray band, or for wave- 
lengths that are large compared to the electron density 
of the scattering medium. For each element, the imag- 
inary and real part of / = /' + if" are related to 
each other through the Kramer s-Kronig relations (e.g. 
IHenke. Gullikson. fc Da"vis"lll993[ D03): 



/o ^ 

w/i(o;) 
2im q r e c 



(3) 



(4) 



where fJ,(u>) is the absorption coefficient at the incident 
energy E — Tiuj. Here the assumption is that for a given 
compound the dielectric function is the sum of the sin- 
gle atoms contributions. However, at the threshold energy, 
the absorption coefficient strongly depends on the chem- 
ical compound. The interaction between a photoelectron 
wave and all the other waves backscattered by the neigh- 
boring atoms creates modulations in the absorption cross 
section. These are in general called X-ray Absorption Fine 
Structure (XAFS). The XAFS features, which are few 
tenth of angstroms wide, have been recently recognized in 
high energy resolution absorption spe ctra of some astro- 
nomical sources fe.g-. lLee et al. ifefJO^ . Only instruments 
with energy resolution better than AA ~ 0.02A can re- 
solve these features. The differential scattering cross sec- 
tion, calculated using the Mie approach, as a function of 
energy is displayed in Fig. ^ for various scattering angles, 
considering the case of Mg2SiC>4 and a grain size of 0.1/mi. 
The presence of fx(oj) in the scattering cross section causes 
spikes at the K-edge of a given element, as seen in the fig- 
ure. We note also that at the softer energies we do not 
expect a dramatic change in the spectral shape of the dif- 
ferential cross section as a function of the scattering angle. 

In our observation, the Cyg X-2 halo is visible down to 
0.4 keV, as the hydrogen column density toward the source 
is relatively low ( Nu ~ 2.17x 10 21 cm" 2 , as measured from 
the Hi emission, iDickev fc Lockman lll 9901 . Therefore we 
need to evaluate the ^ term in Eq.Q]with the Mie theory. 
We then compare the predicted scattered intensity (Eq.HJ, 
calculated using the classical MNR model and the WD01 
model for n(a), with our data. 

3. Data analysis 

Cyg X-2 was obse rved by XMM-Newton EPIC-pn 
I Striider et al~l l200lh in full-frame mode, EPIC-MOS 




Energy (keV) 



Fig. 1. The scattering cross section, calculated using the 
Mie theory, at specific scattering angles vs. energy for 
Mg 2 Si04. A single grain size of 0.1 fim is considered. 



l|Turner et alJl200l|) cameras, and the RGS high resolu- 
tion spectrometers (jden Herder et af~ll200l[) on June 3rd 
2002 for 18.6 ks. The central CCDs of MOS1 and MOS2 
were operated in timing mode and thus contain no imag- 
ing information. We therefore show data of EPIC-pn and 
RGS only. 

We processed the data using the XMM-Newton Science 
Analysis Software (XMMSAS-5.4.1) as well as ad hoc rou- 
tines not included in the standard software. 

3.1. Cyg X-2 EPIC-pn spectral analysis 

A single photon can lead to charges in a single pixel (re- 
ferred to as single events), two neighboring pixels (double 
events), or three or four pixels (if the photon hits a region 
close to the corner of a pixel). Pile- up can occur when 
two photons hit the same pixel in the same read-out cycle 
(energy pile-up) or hit a neighboring pixel (pattern pile- 
up), leading to distortions in the pattern distribution of 
single and double events. The central pixels of the source 
image of Cyg X-2 are heavily piled up. A diagnostic tool 
(XMMSAS task epatplot) applicable to XMM-Newton 
EPIC CCDs makes use of the event pattern distributions, 
which can be precisely modeled as a function of energy 
for X-rays passing through the telescope in the absence 
of pile-up. From a comparison of such a modeled distri- 
bution with a real pattern distribution, the percentage of 
pile-up in a spectrum extracted from a given region can 
be estimated. 

Out-of-Time (OOT) events occur during the read-out of 
a pn-CCD along the read-out direction. As these events 
are accumulated only within a short time and are dis- 
tributed over the whole CCD length pile-up effects are a 
factor of 1000 smaller. After checking with the epatplot 
procedure, we used these events to model the broad band 
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spectrum (0.4-10 keV) of the central source. In princi- 
ple the absolute flux can be evaluated from the OOT 
events. In the case of extreme pile up this is no longer 
possible, due to the "pseudo-MIP" effect: If the charge 
within one pn-CCD pixel exceeds a threshold of about 15 
keV this event is regarded as being due to a minimum 
ionizing particle (MIP) and all events in this CCD col- 
umn and the neighboring columns for t his read-out fram e 
are rejected on board (pseudo-MIPs, iFrevbergl l2003h . 
However, in the case of very strong pile-up this threshold 
can be triggered by normal X-rays. In extreme cases, in 
almost all frames the columns at the center of the PSF 
are rejected. These spatial exposure variations are not 
fully reflected in the event data files and in the XMMSAS 
software. The "pseudo-MIPs" have no appreciable influ- 
ence on the source spectral shape; however, these rejected 
columns affect the determination of the flux measured 
from the OOT events. Indeed, the pseudo-MIP rejection 
occurs preferentially in the columns corresponding to the 
PSF core and therefore in the same columns as the bulk 
of the OOT events, which then get rejected. We modeled 
the spectrum extracted from the OOT events. Note that 
the source position, rather than the recorded position on 
the CCD, was used for the Charge Transfer Inefficiency 
(CTI) correction using XMMSAS. The soft EPIC-pn spec- 
trum is well fitted with a mult i-temperature black body 
for the accretion disk emission l)Mitsuda et alJll984h with 
kT ~0.36keV at the inner radius, plus a comptonized 
black body spect rum for the emission of the neutron star 
( Titarchuk Il994h . The soft (seed) photons have a temper- 
ature fcXb ~ 0.8 keV before being Compton scattered to 
reach a temperature of kT ~ 6.2 keV in an electron cloud 
of thickness r (Tab.^l. We found evidence of an emission 
line at energy E ~ 6.7 keV, consistent with fluor escent 
emission by ionized iron (e.g., iDi Salvo et alJl2002ft . The 
soft spectrum is both absorbed by gas and dust and scat- 
tered by ID, i.e. light is deviated from our line of sight 
"subtracting" photons from the central source spectrum 
(PS95), We implemented in XSPEC a model for the scat- 
tering correction which is based on handy empirical rela- 
tions: r sca = 0.05 x A H - 0.083 and A v = 0.56 A H + 0.23, 
where r sca is the scattering optical depth and Ah is in 
units of 10 21 can . This relation is based on the study of 
25 ROSAT sources (PS95) and it is not critically model 
dependent as long as we are dealing with a relatively small 
correction for r. The optical extinction Ay value is only 
1.3 for Cyg X-2 l|Bradt fc McClintodJlQSafl . therefore the 
influence of scattering in the spectrum is practically neg- 
ligible in the fit for such a low intervening column density 
(PS95). 



3.2. Cyg X-2 RGS spectral analysis 

As the RGS covers only the 0.35-2 keV band, the hard 
component cannot be constrained, therefore the thermal 
comptonization parameters are fixed to the EPIC-pn best- 
fit values. The background spectrum contribution was 



Table 1. Best fit parameters for the EPIC-pn spec- 
trum of Cyg X-2 in the energy band 0.4-10 keV, 
with a disk black body at temp erature kTdb (DISKBB 
in XSPEC, iMitsuda et"aTl Il984h and a compto nized 
spectrum (C0MPTT in XSPE~ iTitarchuk I Il994h . af- 
fected by extinction of gas an d dust (TBABS in XSPEC, 
IWilms. Allen fc McCrav l2000ft . See text for the definition 
of the parameters. Errors are given at 90% confidence level 
for one interesting parameter. 



Nh (xlO^ 
kT db (keV) 
kT (keV) 
kT (keV) 
T 

E Fc (keV) 

(TFe(keV) 

EW Fc (eV) 
X 2 /dof 



0.19 ±0.05 
0.36 ± 0.05 
0.81 ±0.03 
6.2 ±0.1 
1.62 ±0.07 
6.67 ±0.16 
2 +0 - 3 
73 ±53 
1115/1111 



evaluated using the "blank field" observations specific for 
RGS. This guarantees the omission of any halo contam- 
ination in the background. In the case of Cyg X-2, the 
halo contribution to the source spectrum is almost negli- 
gible, since the flux of the source is more than 20 times 
larger than the flux of the diffuse halo. We extracted the 
first and second order from RGS1 and RGS2 for a total 
of four data sets. The data below 7A are affected by low 
and poorly calibrated sensitivity and are rejected. 
Absorption by oxygen in the ISM could be studied in de- 
tail in the spectral region around 0.54 keV. Takei et al. 
(2003) found a complex structure for the oxygen edge re- 
gion in a Chandra-EETG observation of Cyg X-2. In their 
analysis they interpreted the spectrum in terms of ab- 
sorption lines and edges from oxygen in both atomic and 
molecular form. Although the RGS energy resolution is 
approximately 44% less than LETG, we find similar com- 
plexity in the oxygen region: in particular a single oxy- 
gen edge at energy 0.543 keV is an unsatisf actory fit to 
the da ta. We first considered the approach of lTakei et alJ 
(2003J (model 1 in Tab.EJ- We included an additional edge 
in the fit, which improves the fit by A% 2 /A^ = 39/1 (cor- 
responding to a significance higher than 99.5%). The two 
edge energies are fixed: 0.536 (23.13) and 0.543 (22.83) 
keV (A) , corresponding to compound and atomic oxygen, 
respectively. On the other hand, a third edge, strongly re- 
quired by LETG data at 0.549 keV (22.58A, atomic oxy- 
gen), improves our fit only by A% 2 /Af = 5/1 (significance 
97.5%). At an energy of 0.524±0.003 keV, consistent with 
the atomic oxygen Is — 2p transition, an absorption line 
of equivalent width (EW) 1.45j:°; 22 eV is clearly detected. 
Finally, at 0.530 ± 0.003 keV some absorption line-like 
residuals, in addition to the k nown instrumental absorp- 
tion line jde Vries et al. 112003ft . still remain. Including at 
this position a second absorption line in the fit yields an 
energy that is interpreted as the Is — 2p transition of com- 
pound oxygen (Fig.0 upper panel). The measured EW is 
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1.41 ± 0.26 eV, consistent with the findings of lTakei et alJ 
ll2003h . 

However, large uncertainties still remain in the laboratory 
measurements of oxygen bound with other elements and 
the identification of such features i s not conclusive. Other 
laboratory measurements (e.g.. | McLaughlin fc Kirbv I 
Il998t lOorczvca fc McLaughlin 1 1200(1 of atomic oxygen 
around the K edge region would interpret the absorption 
structures as absorption lines from neutral and ionized 
oxygen. This is called model 2 in Tab. The only two 
features in common with the Takei et al. interpretation 
are: the absorption edge at 0.543 keV (22.83A) and the 
Is - 2p transition line at 0.524±0.003 keV. The region 
between these two "standard" features is fit by an ab- 
sorption line at 22.89A, consistent with the ls-3p transi- 
tion of neutral atomic oxygen, and with another absorp- 
tion line which would be consistent with a blend of un- 
resolved lines of Om (at 23.05 A). Finally, the evident 
absorption line at 23.35A, also found by Takei et al., is 
in terpreted as ionized atomic oxyg en (On), as predicted 
bv lGorczvca & McLaughlin"! l|2000l) measurements (Fig. [21 
lower panel). This interpretation (model 2) was also ap- 
plied to Ch andra-'HETG data of a sample of brig ht galac- 
tic sources ijjuett. Schulz. fc Chakrabartv 12004^ . Such an 
ionized component would be interpreted as ionization of 
the ISM, localized in the vicinity of the source. In Tab. 
the results of the two models are shown; there is no sig- 
nificant difference in terms of goodness of fit. 

Table 2. RGS fitting results for the oxygen region in 
Cyg X-2. The parameters for the continuum emission are 
taken from the broad-band spectrum. The energies are 
measured in keV and the Equivalent Width (EW) in eV. 
Model 1 fits the O vicinity with 3 edges (two of which are 
from atomic O) and 2 absorption lines (one from atomic 
neutral O and the other from molecular O). Model 2 in- 
terprets the spectrum in terms of absorption by atomic O, 
either neutral or mildly ionized. Errors are given at 90% 
confidence for one interesting parameter. 
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Fig. 2. Comparison between model 1 (top), and model 2 
(bottom) used to fit the oxygen region, using the RGS 
data. Model 1 includes three absorption edges and two 
absorption lines. Model 2 includes one absorption edge 
and four absorption lines. 

The total equivalent hydrogen column density 
7V H = (2.20 ±0.02) x 10 21 cm" 2 , as measured by the RGS, 
predicts K-shell absorption edges of nitrogen, oxygen, 
neon, and also iron L-shell, as shown in Fig. [3J The 
statistics of the present data allow the determination 
of the physical parameters of the absorption. For each 
element j, we calculated the column density Nj = Tj/o~j 
using the photoelec tric cross section compute d from 
tabulated values llHenke. Gullikson. fc Davis I Il993t 
IVerner fc Yakovlevl Il995(l . The optical depth Tj was 
measured from an absorption edge model. The edge 
model was applied to our best fit continuum, but 
fixing to zero the abundance of the element j in the 
absorption model. Using the IS M abundances listed by 
IWilms. Allen fc McCravl l)2000|) . we derived the equiv- 
alent hydrogen column density. In general we find an 
agreement within the errors between the column densities 
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so derived and the Ah found in our best-fit model 
(Tab. 0). We note that iron shows an overabundance of 
- 20% compared to the lWilms. Allen fc McCravl l|2000|) 
ISM value. 



Table 3. Relevant absorption edges in the RGS spec- 
trum of Cyg X-2. The energy, the corresponding wave- 
length, and the optical depth (r) were obtained from the 
data. The equivalent total hydrogen column density A^h 
was derived from the ISM abundances of Wilms, Allen & 
McCray (2000). These are to be compared with the best- 
fit iV H = (2.20 ± 0.02) x 10 21 cm" 2 , measured by RGS. 
Errors are given at 90% confidence level. 



Element Energy Wavelength 
keV A 



n 



N K-shell 0.409 (fixed) 30.31 0.15 ±0.06 2.7 ±1.1 
O K-shell 0.543 (fixed) 22.83 0.81±0.10 2.8 ± 0.8 
FeL-shell 0.706 ± 0.01 17.56 ± 0.02 0.13 ± 0.01 2.8 ± 0.2 



Ne K-shell 0.869 



+0.02 
0.06 



14.26 



+0.10 
-0.03 



0.06 ±0.03 1.8 ±0.9 





Energy (keV) 



Fig. 4. EPIC-pn image of Cyg X-2 in the 0.5-3 keV band. 
The solid line annuli show three of the extraction regions 
specific for the halo spectral analysis (Fig.0). The circles 
are drawn at 3.3, 4.3, 5.3, 6.3 arcmin. The dashed line 
annulus is the background. These regions were cut by the 
green rectangle across the detector (see § 14.2(1 . For the halo 
spatial analysis fS l4.1[) the annuli were instead logarithmi- 
cally spaced. The barred regions were excluded from the 
data analysis: pile-up region (circle in the center), OOT 
events (box), and the serendipity sources (points). 



relies on (i) an accurate understanding of the instrumental 
scattering by the surface of the gold mirrors of the tele- 
scope and (ii) a careful handling of the pile-up, affecting 
the central part of the source. Pile-up distorts both the 
spectral shape and the surface brightness profile (SBP). 



Fig. 3. ISM absorption edges predicted for ATh~ 2.2 x 
10 21 cm -2 using the best EPIC-pn best fit model, eval- 
uated from the OOT events data. Mg (1.3 keV) and Si 
(1.84 keV) are only marginally evident in the model. 



4. The analysis of the scattering halo 

Thanks to the large effective area of the XMM-Newton 
telescopes, coupled with the large field of view and spec- 
tral resolution of EPIC-pn, the "pure" scattered radiation 
can be extracted from the halo spectrum. Thus for the 
first time we were able to analyze both the spatial and 
the spectral distribution of the scattering halo. The dif- 
fuse emission brightness is a few percent of the central 
source flux. The subtraction of the Cyg X-2 contribution 



4.1. Halo spatial analysis 

As a preliminary step, all spurious sources were removed 
from the data set. The Simbad 1 search reports 60 objects 
in the field of view of Cyg X-2 of which only 3 are known 
X-ray sources. These were not detected in our observation, 
but nevertheless the data at those positions were discarded 
(see Fig.rjJ. The SBP of the scattering halo was evaluated 
extracting the photons at selected energies between 0.4 
and 2keV in bins of 0.25 keV from logarithmically spaced 
annuli centered on the source, excluding the emission com- 
ing from the read-out streak due to OOT events. For the 
background we took an annulus from 9.5' to 13', avoid- 
ing the read-out streak region. The chosen background is 
consistent with the scaled value of the background taken 
from the "blank field" data for EPIC-pn l|Read fc Ponmar] 
l2003h . For the extraction we took into account the diffcr- 

1 http://simbad.u-strasbg.fr 
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Fig. 5. Ratio between single events and all events (up- 
per dash-dotted curve) and ratio between double and all 
events (lower solid curve), for EPIC-pn, around 0.75 keV. 
The dotted lines represent the values expected in absence 
of pile-up. 



Fig. 6. Comparison between the PSF-model (solid line), 
Mrk 421 (asterisks) between 0.9 and 1 keV. For radii 
smaller than ~30", the flatness of the profile is mainly 
due to the mask used in the observation. The counts in 
this region are non zero due to a slight misalignment of 
the mask with respect to the source coordinates in this 
particular observation. 



ent exposures of each quadrant of the detector selecting 
just the time intervals when all four quadrants were on at 
the same time. Any excluded region (serendipity sources, 
OOT events, regions outside the chip boundaries) was of 
course taken into account in the computation of the ex- 
traction area. 

Pile-up is close to 100% at the center of the source. This 
causes the characteristic "hole" in the spatial profile, but 
to a lesser degree also distorts the profile shape up to 
many arcsecs from the source. In order to evaluate the 
dependence of pile-up as a function of the distance from 
the source, we extracted the radial profile of the source 
at different energies for single events, double events and 
the total events, selecting different pattern from the data. 
Where the count rate is low (i.e. no pile-up) the ratio of 
the radial profiles extracted with these different patterns 
(single, double and all events), should be a constant value 
(Fig. 03). We see that for single events pile- up affects the 
profile up to ^40"; thus we studied the data only outside 
this radius. We divided the photon histogram by the ex- 
posure map, and the areas of the annuli, that was also 
corrected for the zones excluded in the photon extraction. 
Moreover, each photon is vignetting corrected by the ra- 
tio of the effective area at the aim point of the telescope 
and at the position where the photon itself is detected. 
The resulting radial intensity distribution is now in units 
of cts/s/arcsec 2 . The next step is the subtraction of the 
PSF from the data. The model for the XMM-Newton PSF 
as a function of energy and off-axis angle was derived by 



the analysis of 110 point-like sources llGhizzardill2002i) . 
The instrumental PSF is described by a King profile: 

PSF oc ( , -. ... ) , (5) 

where r c is the radius of the PSF core, and a the slope of 
the profile. The energy dependence 2 for these parameters 
is a = 1.525 - 0.015 x E and r c = 6.636 - 0.305 x E For a 
full description of the PSF modeling see lGhizzardill|2002j) . 
For each energy we considered in the halo spatial analysis, 
we used the appropriate King profile. We compared this 
PSF profile with the calibration target Mrk 421 (Fig. 0), 
which will be extensively used in this paper as a halo-free 
source, in comparison with Cyg X-2. This bright point- 
like source was observed in a special "masked" mode. In 
this submode the central 11 x 11 pixels (of size 4.128 arc- 
sec each, i.e. a square of ~ 45 arcsec or a half diagonal of 
^0.53 arcmin) were set as bad pixels and thus masked on 
board. This prevents the creation of pseudo-MIPs (§ I3.1fl 
and thus the OOT events can be used to recover the source 
flux value as well as the spectral shape. Thus the Mrk 421 
profile, also in units of cts/s/arcsec 2 , could be easily nor- 
malized to the raw count rate (cts/s) of the source in a 
given energy bin. 

Such a normalization was not straightforward for the 
Cyg X-2 profile, since the central pixels of the source im- 
age are unusable due to the strong pile up. However, since 
the absolute level of the Mrk 421 surface brightness profile 
(in units of arcsec^ 1 ) is known, we could safely normalize 

2 The off-axis dependence of the PSF is neglected as the 
sources we study are located on-axis. 
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the Cyg X-2 halo radial profile. We made the two profiles 
overlap in the interval between ~40"-100", where Cyg X- 
2 does not show significant scattered extended emission 
and therefore the profile is dominated by the instrumen- 
tal PSF. We have seen above that the PSF profile slowly 
changes with energy, hence overlapping the source profile 
with the PSF considering a too large energy band would 
introduce additional uncertainty. For the largest energy 
bin we considered (0.25 keV, § 15 .1(1 the net uncertainty on 
the slope a is ±2%. As the halo analysis is confined to 
radii > 40", the uncertainty introduced by r c (E) is negli- 
gible. We note that in principle, to evaluate the spectral 
flux we could use either the OOT events or the RGS high 
resolution spectrum. The uncertainties in the cross cali- 
bration between PN and RGS flux may reach 20% and 
make the RGS data inappropriate for this purpose. The 
OOT events were also unusable for the normalization due 
to the pseudo-MIP effect (§EHJ- 

At the end of the procedure just described, the resulting 
SBP of Cyg X-2 , which is the summed contribution of an 
extended emission and the instrumental profile, could be 
subtracted by the PSF and then modeled (§ 15. 

4.2. Halo spectral analysis 

For the study of the spatial variation of the spectrum we 
first selected radial distances from the source from 1.3' to 
7.3' divided in annuli of 1'. The background was extracted 
from the same region as for the halo spatial analysis. A 
spatial selection 3 , was applied to ensure that the energy 
response stays constant across the detector. In order to 
have a model independent estimate of how the halo spec- 
tra may change as a function of the distance from the 
source, we normalized the spectra in the rings (in units 
of raw detector counts) to the source spectrum estimated 
from the OOT events fS I3.1(l . In this way, first, all the fea- 
tures belonging to the X-ray emitter itself cancel out and, 
second, assuming that absorption by ID is rather uniform 
within the few arcmin across the halo, also the absorption 
component is eliminated from the halo spectrum. The re- 
sult of this procedure is shown in Fig. The spectra of 
the Cyg X-2 halo are plotted in order to illustrate the be- 
havior of the halo with increasing angular distance from 
the source, as the PSF contribution becomes less impor- 
tant. The extraction regions of these spectra are marked in 
Fig. Here the inner radii range from 3.3' to 5.3', where 
the halo is more relevant. Note that at this stage, the in- 
strumental mirror scattering spectral energy distribution 
is still to be subtracted. The vertical axis of Fig. is in 
arbitrary units as the absolute normalization of the source 
spectrum cannot be recovered from the OOT events (see 
§ 13.1(1 . Below 2 keV, a decline of the curve is evident and 



< 0.001 - 




3 the RAW CCD coordinate selection was: RAWY=180 - 
199. This RAWY range always contains the source if it is ob- 
served on-axis and is well calibrated, and allows us to avoid the 
uncertainties related to the RAWY dependence of the energy 
response. 



Energy (keV) 



Fig. 7. The ratio between the halo spectra and the central 
source spectrum. At this stage the data are still convolved 
with the PSF. The spectra were extracted from annuli - 
from top to bottom- centered at 3.8, 4.8, and 5.8 arcmin 
respectively (Fig.0J). At larger radii the PSF contribution 
becomes less important. Dashed vertical lines: energies of 
relevant ID elements (O, Fe, Mg, Si). Error bars are 68% 
confidence level (la). The vertical axis units are arbitrary 
as the absolute flux of the source could not be recovered 
from the OOT events (§E3J- 

the presence of three absorption features corresponding 
to the energy of O (0.54 keV), Mg (1.3 keV) and Si (1.84 
keV) are observed. The vertical dashed lines indicate the 
position of these features and for completeness, the Fe L 
edge position (0.709 keV) is also marked. In the absence of 
any extended emission (i.e. for a point-like source) the re- 
sult of the above exercise on the spectrum extraction is the 
pure instrumental mirror scattering. This was done for the 
point source Mrk 421, to which exactly the same proce- 
dure was applied (Fig.[S]for an example at radius 180"). 
Since the point-source data are real data (though with 
high statistics), they are still affected by noise. Thus, to 
avoid loss of information on the halo spectra, the Mrk 421 
spectrum was smoothed before subtracting it from the 
data. The PSF spectrum does not change significantly as a 
function of the scattering angle at soft energies, keeping its 
almost constant shape between 0.4— 1.2keV. As we do not 
have precise information on the absolute flux of Cyg X- 
2, in order to perform the PSF subtraction we made the 
halo spectrum collected at each annulus coincide with the 
homologous spectrum of Mrk 421, in the spectral region 
where the emission is due to scattering by the telescope 
mirrors only (> 3.5 keV). The PSF-subtracted scattered 
radiation is shown in Fig.^Jfor an observed angle of 4.8'. 
From the lowest observed energy up to the silicon edge 
energy, the PSF subtraction is not a concern, as the ex- 
tended emission is well above the PSF spectrum for a wide 
range of scattering angles (see Fig. El for an example at 1 
keV). The comparison PSF- halo is less straightforward for 
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Fig. 8. XMM EPIC-pn pure-scattering mirror spectrum 
at a radius 180" from the source, extracted from the halo- 
free source Mrk 421. We note a smooth increase of the 
scattering above 2 keV for larger radii. 



the Si feature and the continuum just below 1.84 keV, as 
the scattered emission becomes weaker with respect to the 
PSF. In this region, a small difference in the chosen nor- 
malization may affect the results of the halo modeling. 
The PSF subtraction has the effect of making the Si fea- 
ture appear even deeper than what already observed from 
the raw data (Fig. 0. If, for example, the normalization 
of the halo spectrum were ±10% of what we calculated, 
the error bars for the Si feature in the PSF-subtracted 
spectrum (Fig. Ill|) would increase of an additional 6%, 
considering that the ratio between the PSF and the raw 
data at this position is ~ 30% at o bs = 4.8'. 

5. The halo modeling 

There are two critical terms in Eq. ^ One is the differen- 
tial scattering cross section and the other is the assumed 
grain size distribution. For both the spatial and the spec- 
tral modeling of the scattering halo, we used the differ- 
ential scattering cross section calculated using the Mie 
theory. For each value of the parameter X — 2ir a/ A, 
where a is the grain size and A is the wavelength of 
the incoming wave, we calculated the parameter da/dfl 
for 115 scattering angles 9 sca from to 3°. The scat- 
tering cross section depends also on the diffraction in- 
dex m which is determined by the material and is en- 
ergy dependent. In this study we use the tabulated values 
of m for graphite, olivine, and pyro xene as determined 
by^Ienkc. Gullik son. fc Davis I l|l993() except in one case 
(8 lfi.l[l . when we test the specific ID composition proposed 
in D03. In that case, the calculation for m takes into ac- 
count the X AFS near the edg es. 

The code 4 l|Wiscombe I Il08c5f) used to generate the scat- 

4 ftp:/ /climate. gsfc.nasa.gov/pub/wiscombe/Single_Scatt/ 
Homogen_Sphere /Exact JVlie / 



tering cross sections is reliable for X < 2 x 10 4 . For large 
values of a (say a > 1.2/im for E = 2 keV), the anomalou s 
diffraction theory should be used l|van de Hulst I Il957|) . 
We igno re very large grains in our calculation since, as 
shown in lDraine fc Tanl ((2003) , radii > 0.4/^m contribute 
less than 1% at large scattering angles, and less than 20% 
at ^sca < 100". In bright sources, like Cyg X-2, pile-up 
hampers the possibility to investigate the halo at small 
scattering angles, where the effect of very large grains, or 
grains located very clos e to the source may be relevant 
l|Predehl fc Klose Ill99q) . We also restricted our modeling 
to energies <2 keV, above which the halo contribution 
drops dramatically in the case of Cyg X-2. 
Thus, the chosen grain size interval is 0.005-0.25 /im or 
0.00035-0.8 (im when the MRN dust size distribution 
model or WD01 model is adopted, respectively. The dust 
size intervals were divided in 200 logarithmically spaced 
size bins. We allowed the power law index of the MRN 
distribution to vary by 20% around the typical value 3.5. 
WD01 tested their grain size distribution for two dif- 
ferent values of the ratio of the total over selective op- 
tical extinction: i?y=3.1 and 5.3, and for different car- 
bon abundances. Cyg X-2 is located at galactic latitude 
b = -11.3° where the ISM is diffus e (no CO detected, 
iDame. Hartmann. fc Thaddeudl200l|) . thus we considered 
Rv=3.1. We used the set of parameters for slopes and 
coefficients of the dust distribution corresponding to a 
carbon abundance in PAH alone of 6xl0 -5 (Tab. 1 of 
IWeingartner fc Draine ll200l|) . For both models, the lower 
and upper limit of the integral on the dust distribution 
parameter x were left as free parameters. 

5.1. Spatial Modeling of the Halo 

The SBP was extracted and subtracted from the PSF con- 
tribution as described in S I4.1I At a fixed energy, the model 
has three free parameters (Eq.^l: da/dfl, x, and n(a). The 
best fit was reached through x 2 minimization. We con- 
sidered rays only scattered once before being observed. 
Double scattering occurs for optical depths r sca close to 
unity, indicating a very high dust column density (PS95, 
ICostantini fc Predehl I [2005) , which is not observed to- 
ward Cyg X-2. In Fig. M the Cyg X-2 SBP of the halo 
at 1 keV is shown. We tested the MRN and the WD01 for 
the dust size distribution. Both models provide an accept- 
able fit in terms of % 2 (x^ ed =1.29 and 1.33, respectively) 
. The WD01 distribution spans a wider range of grain 
sizes. In particular, scattering by grains with size a in the 
range 0.25 — 0.4/im have the effect of enhancing the halo 
at smaller radii ( ^ 200"). The intensity of the halo is pa- 
rameterized by the scattering optical depth r sca , defined 
as: Zf rac = Jhalo/itot = 1- e~ Tsca , where /halo is the flux of 
the scattered emission, and /tot is the total source emis- 
sion (PS95). At 1 keV we measured r sca = 0.054±0.018 
and 0.067±0.018 for the MRN or the WD01 model, re- 
spectively. The error quoted here is statistical and does 
not include any uncertainty in the background subtrac- 
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tion. The best fit shows a minimum and maximum value 
for the fractional path at which the scattering occurs, x, 
of ~0.001 and ~0.6, respectively. In this interval the halo 
profile shows a smooth distribution of dust along the line 
of sight up to a fractional distance x ~ 0.6, after which 
the halo is unaccessible due to pile-up. However, if the up- 
per limit of x is constrained to be close to 1 (we fixed it 
at 0.99) indicating that we are actually observing scatter- 
ing occurring at all distances, the fit worsens for both the 
WD01 and MRN models (Ax 2 = 37 for Au = 1, corre- 
sponding to a significance > 99.5% ). We then extended 
this analysis to the energy range at which the halo is ob- 
servable. In Fig. ^| we show the total scattering optical 
depth, derived from the SBP using energy intervals of 0.25 
keV, in the energy range 0.4-1.9 keV. The large bin size 
smoothes out any features, leaving just the general shape 
the spectral energy distribution of the scattering optical 
depth. The solid line in Fig. ^| refers to the theoretical 
value of r sca as predicted by D03 at the mean energy of 
the extraction bin. This was derived by multiplying the 
theoretical value of cr sca (D03) with the hydrogen column 
density toward Cyg X-2 that we measure. 
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Fig. 10. The scattering optical depth as a function of 
energy, as measured from the SBP of Cyg X-2 halo for the 
MRN model (filled squares) and WD01 (empty squares). 
The data were collected in energy bins 0.25 keV wide. The 
solid line refers to the value of the total r at the center of 
the bin, as predicted in D03. 




Radius (arcscc) 

Fig. 9. The data (halo+PSF) (asterisks) compared to the 
PSF (dotted line) around 1 keV. The dashed and the solid 
thin lines are two different halo models (MRN and W01), 
while the solid and dashed thick lines indicate the best- 
fit to the total data (model+PSF) relative to MRN and 
WD01, respectively. 



5.2. Spectral Modeling of the Halo 

In Fig.^Jwe show the pure scattered radiation energy dis- 
tribution (i.e. PSF subtracted), collected in a ring centered 
at 4.8'. The subtraction of the mirror scattering and of the 
central source contribution was described in § 0] We avoid 
the data below 0.4 keV as the effective area calibration 
may still suffer from uncertainties in this range. At soft en- 



ergies, the halo spectra collected in different annuli, shown 
in Fig. \7\ do not vary dramatically (within statistical er- 
rors) among each other after the PSF subtraction. Above 
2 keV the differential scattering cross section does change 
dramatically as a function of the scattering angle (Fig. ^) ; 
however the contribution of the PSF wings is dominant at 
this energy and this change is unobservable in the present 
data. In principle, after the PSF subtraction, the model- 
ing of the halo spectra, at each angular distance, would 
reveal the local properties of the dust grains. Indeed, in- 
homogeneity in the dust spatial distribution, size of the 
dust particles, and chemical composition would result in 
a spectral change. Unfortunately, these changes are tiny 
with respect to the instrumental uncertainties. The model- 
ing here refers just to the ring around 4.8', where the halo 
is brighter and thus the statistics are highest. The spec- 
tral features are still very well observable after the PSF 
subtraction. If artificially fitted with Gaussian profiles, the 
significance is 2.3ct, 3.2ct, and 4ct for O, Mg, and Si, respec- 
tively. Although the dust features of oxygen, magnesium 
and silicon are clearly detected in the scattered emission, 
still the statistics do not allow a unique interpretation of 
the data. In the attempt to give a quantitative picture of 
the observational evidence, we chose a model with only the 
dominant well-known compounds of the diffuse interstellar 
dust environment. We calculated the scattered intensity 
expected at the chosen angular distance, for the most com- 
mon constituents of interstellar grains, i.e. graphite (car- 
bon) and silicates. Silicates are found mostly in the form 
of Mg 2;c Fe2(i_ x )Si04 and MSi0 3 , where M is either Mg 
or Fe, thus we included olivine in the three different forms 
of MgFeSiG-4 (p=3.8 g cm~ 3 ), Fe 2 Si0 4 , (p=4.39 g cm" 3 ) 
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and Mg2Si04 (,0=3.27 g cm 3 ). Pyroxene are in the form 
of MgSi0 3 (p =3.2 g cm" 3 ) and FeSi0 3 (p=3.8 g cm" 3 ). 
For a given dust size distribution model (MRN or WD01), 
the relative contributions of the different grains materials 
are left as free parameters. Moreover, we considered also 
a model where only MgFeSi0 4 is taken as representative 
compound for silicates (D03). For this model, the dust size 
is distributed according to WD01, and in the scattering 
cross sections XAFS are included. The goodness of the fit 
was evaluated by minimizing \ 2 - 

The data are ignored below 0.4 keV, therefore we do not 
expect to put any constraint on the carbon component, 
whose main feature is at 0.28 keV. However, the carbon 
component is included in the fit since it may influence 
the shape of the continuum at energies above the car- 
bon edge. Carbon contribution was allowed to vary be- 
tween 20% and 30% of the total amount of dust ijWhittet I 
l2003h . for all the models we tested. The data are mod- 
eled by a mixture of silicate and pyroxene. In Fig. 1111 we 
display the three dust models tested. A mixture of sili- 
cate compounds (labeled MRN and WD01) seems to ac- 
ceptably interpret the data. C, O, Fe, Mg and Si alone 
account for 95% of the dust components. All elements ex- 
cept O account f or 15-30% of th e total amount of dust, 
in various forms llWhittetlEool . Here we assume that 
these elements describe 100% of the observed scattering. 
The linear combination of the compounds contribution is 
shown in Fig. [H for the MRN distribution. In the WD01 
case, the combination of compounds is not significantly 
different from MRN fFig. Ill[). The relative contribution 
to unity is Mg 2 SiO 4 =0.42, C=0.26, FeMgSiO 4 =0.054, and 
FeSiO 3 =0.25. Fe 2 Si0 4 and MgSi0 3 contribute for a negli- 
gible fraction (<0.001). Counting the contributions of the 
single elements, we obtain that roughly 42% of O, 26% of 
C, 11% and 13% of Si and Mg, and 5% of Fe are needed 
to fit the data. 
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Fig. 11. EPIC-pn data of the Cyg X-2 halo, extracted 
at ~4.3' from the source and PSF subtracted (asterisks), 
compared with the best-fit models: WD01 (dashed dot- 
ted line), MRN (solid line), and D03+WD01 model (dia- 
monds). O-K, Mg-K, and Si-K edge energies are indicated 
by dotted lines. 





6. Discussion 

6.1. The chemistry of dust grains 

For the first time the signature of the elements locked 
in dust grains and responsible for the scattering of the 
X-rays is detected. The best fit indicates a major contri- 
bution by olivine and pyroxene. We cannot exclude the 
presence of other different compounds. Being the best fit 
a linear combination of scattered intensities for a given 
compound, adding too many components would not be 
necessarily a true physical interpretation of the data. In 
principle, the depth of the spikes in the scattered spectrum 
tells us about the intrinsic properties of the dust grains. 
Magnesium and silicon are not detected but marginally 
in absorption (Fig. |3J), while they are prominent features 
in the scattered spectrum, even in the raw data at the 
angles where the statistics is maximal (Fig. 0. Only a 
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Fig. 12. Relative contributions to the best-fit model, us- 
ing the MRN distribution. The data are the same as in 

Fig. EH 



deeper observation will allow us to quantify this possible 
discrepancy by comparing the column densities for scat- 
tering and absorption derived for Mg and Si. As Fig. 
shows, for an absorbing equivalent hydrogen column den- 
sity of 2.2xl0 21 cm" 2 , the iron L-shell at 0.706 keV is 
clearly measured. In the scattered spectrum iron is hardly 
detectable. The best fit model ("Fig. 1121) requires a certain 
amount of Fe, although not precisely quantified, locked in 
silicates, with a Mg to Fe ratio of ~5: (1.4-2.5). The deple- 
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tion of Mg, as well as Si, depends on the density of the 
dust environment. Differently, Fe is found to be highly de- 
pleted (80-100%) from the intercloud medium to the dens- 
est environments. Thus, we can suppose that Mg is more 
abundant in silicates than Fe. In particular, a mixture of 
olivine and pyroxene with a ratio of 5:2 for Mg: Fe, would 
totall y account for the depletion of Mg and Si ijWhittet I 
120031) . In this frame, less than half of the available iron 
grains are locked in silicates and the rest in other forms. 
The best fit of the scattered halo at ^4.8', within errors, is 
in agreement with this simple prediction. If the lower limit 
is taken, the re would be more room for iron to be locked in 
other forms. IWeingartner k. Draine I l)l999|) suggested, for 
instance, that up to 60% of Fe could be associated with 
graphite to form very small grains (a ~ 10 — 15 A). A 
Mg:Fe ratio of 1:1, as prescribed by the D03 dust mixture, 
fails to interpret the spectral energy distribution of the 
dust toward Cyg X-2. This again reinforces the idea that 
more compounds containing Mg and Si should play a role 
in the scattering. As shown in Fig. 1121 the data requires 
a major contribution of Mg2SiC>4 rather than iron com- 
pounds like FeSiOa . This dust mixture leads to a stronger 
contribution of Mg with respect to the D03 model. 
This interpretation is influenced by the PSF subtraction, 
which can artificially enhance the depth of the Si fea- 
ture (but not significantly for Mg which lies at 1.3 keV, 
where the halo is well above the the PSF) and the con- 
tinuum around it (§ I4.2fl . Another unknown is the role 
of carbon. If the contribution of carbon between 0.4 and 
2keV varies substantially from the 20—30% of the total 
budget, then the contribution of silicate would need to 
be revisited. The quality of spectral data does not al- 
low to appreciate a substantial discrepancy between the 
WD01 and MRN dust size distribution. As discussed be- 
low I6.3fl . at the scattering angles we consider for the 
spectral analysis (around 4.8'), the bulk of the scatter- 
ing is caused by grains with quite "standard" sizes and in 
this case the MRN and WD01 do not differ dramatically. 
However, with a deeper observation we would differenti- 
ate more between the two models at each scattering angle 
and this would influence the relative contribution to the 
scattering of the different compounds. We have seen that 
the modeling of Si suffers from the highest uncertainty 
in this analysis and we cannot draw any conclusion on 
the basis of this feature. However, the prominence of the 
Mg feature suggests a significant contribution of magne- 
sium compounds. The oxygen region of the scattered spec- 
trum is well interpreted by D03. This is because of the 
dielectric functions, which include XAFS near the edge 
energy (D03). When convolved with the spectral resolu- 
tion of these data, such sub-structures are almost totally 
canceled: the difference in depth b etween the D03 and the 
iHenke. Gullikson. fc Davis"! l|l993() dielectric functions are 
indeed ~1% for Mg and Si, and ^20% for iron. Oxygen is 
the only element where the discrepancy is noticeable (~ a 
factor two). Evaluating the contribution of the single ele- 
ments, assuming that C, O, Fe, Mg, and Si are causing the 
totality of the scattering, we find that oxygen can be com- 



pletely explained in terms of silicate. This supports the 
idea that oxygen is preferably locked in these materials, 
leaving little room to other O compounds (e.g. OH, H2O, 
IWhittet etai~ll200lh . at least in the diffuse ISM. The fea- 
tures of the scattered halo, as the one shown in Fig . 1111 
were predicted by the scatte ring theory (e.g. iHavakawal 
ll970tlPredehl fe Klose lll99l D03). The dust physical pa- 
rameters extracted from our modeling seems to bolster 
this interpretation. Another possibility for producing a 
complex halo spectrum would be an uneven absorption 
at an angular scale of arcminutes around the source. If, 
for example, the absorption toward the central source is 
few percent less than at the outer parts, the halo/source 
ratio will show extra absorption at soft energies, similar in 
shape to what we observe (Fig.EJ). In this case, at the en- 
ergy of the absorption edges, we would expect also broad 
and asymmetric residuals. However, we found the spectral 
shape of the halo to have circular symmetry around the 
central source and this is in conflict with a clumpy struc- 
ture of the medium, and the features we observe in the 
scattered spectrum seem not to have an asymmetric, edge- 
like, shape. Moreover, there should be regions where, on 
the contrary, iVu is lower than toward the central source. 
In this case, the halo/source ratio should show an excess 
at soft energy and this was not observed. Although the 
assumption that in the diffuse ISM dust is homogeneously 
distributed is surely simplistic, in the case of Cyg X-2 we 
cannot prove that spatial variations of Nh on arcminutes 
scales play a major role in the halo shape. 

6.2. Oxygen in the ISM 

The most prominent absorption features in the RGS band 
arise from oxygen transitions. The Oi resonant line cor- 
responding to the ls-2p t ransition is clearl y detected. 
Whether we can observe O l|Takei et al.ll2003j) in molecu- 
lar form is not obvious. The uncertainty is enhanced con- 
sidering the lack of laboratory measurements for oxygen 
locked in compounds at these energies. The RGS resolu- 
tion does not allow us to give a unique interpretation of the 
spectral region at energies higher than the oxygen edge. 
In terms of x 2 , the addition of one or possibly two more 
absorption edges or, alternatively, modeling the spectrum 
with three absorption lines (Oi ls-3p, On ls-2p, Oin, ls- 
2p), pro vides similar results. From the curve of growth 
analysis, I.Tuett. Schulz. fc Chakrabartvl 1)200411 calculate 
that the contribution of the ionized oxygen should be 10% 
of the neutral phase. The mildly ionized oxygen would 
then be produced by charge exchange with Oi. However, 
in principle, absorption by oxygen in compounds is ex- 
pected. In the present XMM-Newton observation, in ad- 
dition to absorption, we could also study the scattering 
by ID and this allows us to state that oxygen locked in 
dust is actually present on the line of sight of Cyg X-2 
(Fig. II 1|) . The modeling of the halo spectrum shows that, 
to first order approximation, the oxygen feature is well 
explained if all the oxygen available in dust is locked in 
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silicates. This transla tes into an oxygen to hydrogen ratio 
O dust /H= 180 ppm l|Cardelli et al. Ill996|) in the line of 
sight of Cyg X- 2. If the oxygen abundance O/ H is taken 
to be 490 ppm l|Wilms. Allen fc McCra71l2000h . then the 
depletion value for O, defined as the gas to total ISM ratio, 
is 0.63. With the total hydrogen column density derived 
from the RGS measurement (Nh ~ 2.2 x 10 21 cm -2 ), we 
can then predict an absorption column density for oxy- 
gen locked in (silicate) d ust. We obtain N ^„ for dust 
- 3.85 x 10 17 cm- 2 . In the lTakei et"ai] l|2003h model (col- 
umn 1, Tab.[5J), the edge attributed to oxygen in dust has 
an absorption optical depth of 0.27 ± 0.06 at 0.536 keV. 
Using the absorption cross section aSL = 1-78 x 10~ 19 cm 2 , 
(D03, iLi et al. Ill995h . the column density for this edge is 
1-51q'4 x 10 17 cm -2 , which is too low compared to the the- 
oretical value we obtained above. One possibility is that 
other absorption edges by oxygen in solid form are missing 
the detection. Indeed the energy and the structure of the 
edges slightly changes depending on which compound of 
oxygen and silico n is considered fe.g. lOwens et al.l l2002: 
IWong et al .111 9951) . and therefore some of them may be 
blurred and difficult to detect with the present resolution. 
This calculation strongly depends on the adopted value of 
°"abs- Here we considered the peak abso rption cross sec - 
tion for oxygen locked in Mg 2 Si0 4 (D03, ILi et al. lll995h . 
as determined by D03, with the underlying assumption 
that the oxygen edge profile is the same as the Si pro- 
file. However, small variations of <r® hs , lead to substan- 
tial variations of the value of Nh- In fact, if the value 



Indeed these two conditions have the same effect 
of producing very narrow halo components in the profile. 
The data modeling of both the halo spatial profile and 
the halo spectral distribution show that both the MRN 
and WDOl dust size distribution can be applied despite 
the different amount of dust predicted for different scat- 
tering angles (Fig. EJ). Grains with sizes > 0.25/im play a 
more important role in the WDOl model for radii ^ 200" 
at 1 keV. In the size range 0.001-0.25 fim, the MRN 
and WDOl mo del do not differ dramatically for i?y=3.1 
(see Fig. 2 of IWeingartner fc Draine ll2001^ and for the 
grain sizes roughly between 0.1-0.2/im typically produce 
the bulk of the scatt ering halo at angular radii between 
100" and 1000" (e.g, IMathis fc Lee1ll99lHDraine fc Tanl 
2003), which is the region that could be directly studied 
in the SBP of Cyg X-2. The contribution of small grains 
(a 0.05) would be best investigated if we could access 
the region beyond 1000", but unfortunately, in this ob- 
servation the halo begins to fade, making the modeling 
challenging beyond ^650". 

The value of the scattering optical depth derived from 
the SBP at 1 keV (r sca = 0.054 ± 0.018) is larger, but 
consistent within the errors, with the value derived by 
the ROSAT halo (r sca ~ 0.039, PS95) obtained at -1.06 
keV, if the MRN distribution is used. The PS95 model- 
ing indeed started also from the MRN model, but leaving 
free some parameters in it, among which the maximum 
grain size and the slope of n(a). Their best fit requires 
a maximum grain size of only 0.15 /an. This makes the 
^abs i s c l° ser t° what predicted for atomic oxygen (a = halo model flatter and the derived value of r sca lower than 

what we measure. The WDOl distribution predicts instead 
almost the double (r sca = 0.067 ± 0.018) of the ROSAT 
result. This trend is visible at all energies fFig. HOfl . This is 
mostly due to the increased scattering "power" by larger 
grains at smaller radii predicted by the WDOl (Fig. |§J. 
When the total scattering optical depth (for both MRN 
and WDOl case) as a function of energy is compared 
with a theoretical model (D03, Fig. I10|) . we see that it 
does not measure a substantial part of the scattered ra- 
diation. Indeed the model over-predicts the data at all 
energies apart from perhaps the last point. The distri- 
bution of r sca flattens toward lower energies, while the 
model, is significantly steeper. As noted in D03, the dis- 
crepan cy was also found for other halo analysis (e.£ 
PS9 5. ISmith. Edgar, fc Shafer I l2002t IWoo et alJll994 
The IPraine fc TanT rchfW measurement instead has been 
found in agreement with the model (D03). 
The value of Jf rac (and thus of r sca ), extracted at each 
energy, is model dependent and, moreover, can be sub- 
stantially influenced by instrumental effects. For exam- 
ple, the scattering angle range accessible to this study 
(200"— 600") strongly privilege the o bservation of the scat- 
tering of 1 - 2 keV photons (e.g., IMathis fc Leelll99ll 
D03), while the bulk of the emission for softer photons 
peaks at larger scattering angles. The estimated values 
for the total scattering optical depth, and in particular 
the ones related to the soft scattered photons are certainly 
lower limits. A significant fraction of the halo, in the form 
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4.8 x 10~ 19 , Verner & Yakovlev 1995; a = 3.85 x 10" 
Henke, Gullikson & Davis 1993), then the column density 
associated to the 0.536 keV edge reconciles to the value 
we obtained from the scattered radiation analysis. As a 
cautionary note, the effective area of EPIC-pn diminishes 
toward lower energy (60% less at 0.54 keV than at IkeV) 
and the statistics at the oxygen feature may prevent us 
from a full description of the dust contribution to the to- 
tal extinction. Moreover, the value for the total O /H ratio 
is not yet solidly established ijjensen et al. 2005, and ref- 
erences therein). Finally, the knowledge of the theoretical 
wavelengths of these dust absorption features is very frag- 
mentary, making their interpretation somehow aleatory. 
Taking the N^ bs derived starting from the scattering anal- 
ysis (and t = N^gCr^g = 0.68, if D03 value is taken 
for a) as an upper limit, absorption by oxygen in dust 
should be visible in absorption (edges and lines) using the 
RG S, and instrum ents with even higher energy resolution 
l|Takei et al.ll2003[h 



6.3. The dust size distribution 

In this XMM-Newton observation, the pile-up prevents 
us from studying the halo at angular radii smaller than 
~40". This observational limit also hampers any detection 
of scattering either by large grains or s tandard grains lo- 
cated very close to the emitting source l|Predehl fc Klose 1 
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of fainter or narrower components, may be masked by the 
PSF wings (at small radii) or unaccessible because of the 
faint scattered emission with respect to the background 
(at larger radii). 

6.4. The dust distribution along the line of sight 

The distribution of dust along the line of sight could be 
studied through the SBP of the halo. The distribution of 
dust seems to be evenly distributed, at least for a frac- 
tional distance of the total path x < 0.6, corresponding 
to a linear distance between 4.3-6.7 k pc, (depending on 
the distance estim ates for Cyg X-2; lOrosz fe Kuulkersl 
Il999t ISmald IT998|h However, if the dust distribution is 
imposed to be uniform up to a fractional distance of 0.99 
the fit worsens significantly suggesting that our line of 
sight passes through different dust clumps. A contribution 
from scattering events closer to the source is indeed likely 
to be present, but fainter or na rrower halo components, 
as detected in other sources (e.g.lSmith. Edgar, fe Shafer I 
Eool lDra,ine & Tan 11200.4 ICostantini fc PreHehl IbOO.fl . 
cannot be investigated using XMM-Newton for radii < 12" 
(and in any case not for Cyg X-2 because of the pile-up). 

7. Conclusions 

We have presented XMM-Newton results on the effect of 
scattering and absorption by ID along the line of sight to 
the bright X-ray binary Cyg X-2. This study led to the 
unprecedented detection of the elements in the ID respon- 
sible for the X-ray scattering: oxygen, magnesium, and 
silicon. To first order, the modeling of the pure scattered 
radiation suggests a major contribution of silicates in the 
form of olivine and pyroxene, in the energy range 0.4-2 
keV. The best fit of the scattered spectrum shows that 
the ratio of Mg to Fe, locked in dust grains, is 5:(1.4:2.5). 
This is consistent with a picture in which Mg and Si are 
for the most part locked in silicates IjWhittet I 2003) . 
The contribution of carbon, a fundamental constituent of 
ID, could not be quantified as its most prominent feature 
(0.28 keV) lies below our selected EPIC-pn energy band. 
In the RGS spectrum, we studied the complexity around 
the oxygen edge, investigating the possibility of absorption 
by atomic and molecular oxygen, as suggested by Takei et 
al (2002), in comparison with absorption by atomic neu- 
tral and ionized oxygen ij.Tuett. Schulz. fc Chakrabartvl 
2004|. The RGS resolution is not sufficient to prefer one 
interpretation to the other, but the complementary infor- 
mation from the EPIC-pn analysis of the scattered spec- 
trum allows us to detect oxygen locked in dust, preferably 
in the form of silicates. Starting from the scattered halo 
spectral modeling, we estimated that the absorption col- 
umn density we expect by oxygen locked in silicates is 
indeed measurable using the RGS. The value of this col- 
umn density seems too high compared to what measured 
for the absorption edge that iTakei et al.l l|2003h interpret 
as arising from dust. However, instrumental/theoretical 
uncertainties makes this result not conclusive. 



The study of the spectral energy distribution of the scat- 
tered radiation, performed in the halo region where the 
signal-to-noise ratio was best (around ~4.8'), stressed the 
need of using an accurate theoretical approach to the 
data. The full Mie theory had to be used to model the 
data satisfactorily, especially below 2 keV where the chem- 
istry of the halo can now be studi ed. This approach was 
already applied to R OSAT data (|Smith fc Dwek I ll99St 
iDraine fc Taii1l2003h . With XMM-Newton we could ex- 
tend this analysis, performing spatially resolved spec- 
troscopy of the halo that could not be interpreted unless 
the Mie differential scattering cross section was used. 
The modeling of the SBP shows that the dust is uniformly 
distributed along the line of sight at least for a fractional 
distance of the total path x < 0.6, corresponding to a lin- 
ear distance between 4.3-6.7 kpc, depending on the source 
distance estimates. However, a uniform dust distribution 
along the complete path toward Cyg X-2 is not required 
by the data, hinting to a clumped structure of the dust 
for x > 0.6. Within the instrumental uncertainties, the 
data are acceptably fit by both a MRN and a WD01 dust 
size distribution. The inferred scattering optical depth is 
approximately 0.054 and 0.067 at 1 keV for the MRN and 
WD01 distribution, respectively. We extended the mod- 
eling of the SBP of the halo to the 0.4-2 keV band. The 
derived values of the total r sca as a function of energy are 
systematically lower than what predicted by the theory 
(D03), pointing out that some halo components may be 
easily missed due to instrumental effects. 
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